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Abstract 



We present quantum (and classical) Bianchi I model, with free massless scalar field, 
of the Universe. Our model may be treated as the simplest prototype of the quantum 
BKL (Belinskii-Khalatnikov-Lifshitz) scenario. The quantization is done by making 
use of the nonstandard Loop Quantum Cosmology (LQC). Since the method is quite 
new, we present in details its motivation and the formalism. To make the nonstandard 
method easily understandable, we include its application to the FRW model. We solve 
the Hamiltonian constraint of the theory at the classical level and find elementary 
observables. Physical compound observables are defined in terms of elementary ones. 
We find that classical Big Bang singularity is replaced by quantum Big Bounce transi- 
tion due to modification of classical theory by holonomy around a loop with finite size. 
The energy density of matter fields at the Big Bounce depends on a free parameter A, 
which value is expected to be determined from future cosmological observations. The 
phase space is divided into two distinct regions: Kasner-like and Kasner-unlike. We 
use the elementary observables to quantize volume and directional volume operators 
in both cases. Spectra of these operators are bounded from below and discrete, and 
depend on A. The discreteness may imply a foamy structure of spacetime at semi- 
classical level. At the quantum level an evolution of the model is generated by the 
so-called true Hamiltonian. This enables introducing a time parameter valued in the 
set of all real numbers. 
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Introduction 



It results from cosmological observations that the Universe emerged from a state char- 
acterized by extremely high density of matter fields called Big Bang or cosmological 
singularity [U [21 Gl E] ■ The observational data are in comparatively good agreement 
with calculations obtained within the standard cosmological models (SCM). The lat- 
ter is based on General Relativity (GR), and includes dark matter, dark energy and 
the inflation transition. However, SCM provides a phenomenological description. 
Deep understanding of the cosmological data, in particular of an early Universe re- 
quires an elementary quantum theory being able to explain creation of spacetime 
filled with matter fields. 

Probably, the best starting point for finding such a theory is quantization of the 
BKL (Belinskii, Khalatnikov, Lifshitz) scenario [3 El [7]. First of all, it is clear from 
this scenario that a general solution of the Einstein equation with the cosmological 
singularity does exist. By this term we mean a singularity in time on spacelike hyper- 
surface characterized by blowing up the curvature invariants together with diverging 
energy densities of matter fields. The BKL solution of GR is general and stable. By 
general we mean containing a non-zero measure subset of all initial conditions. Sta- 
bility means that no infinitesimal perturbations of initial conditions is able to change 
the singular character of the solution. 
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Secondly the BKL scenario says that in asymptotic vicinity of the singularity this 
general solution has complicated oscillatory behavior of chaotic character. There is 
considerable support for this scenario both from analytical [8j [9] and numerical [TO] 
investigations. What is important in the BKL the dynamics at any spatial point can 
be approximated by that of the homogeneous (but in general non-isotropic) models 
which are called Bianchi models [SI El E]- 

In the BKL scenario, time derivatives of gravitational field is shown to dominate 
over spatial derivatives for long stretches of time. During such periods called Kasner 
epochs, an evolution of gravitational field may be approximated by the Bianchi I 
model [TTj . 

In each Kasner epoch the Universe is approximated by the Bianchi I metric with 
some specific set of parameters kj satisfied equations Y^=i k* — 1 an d Y^=i k^+k^, = 1 
where describes density of matter fields. There exist two classes of solutions. The 
first, called Kasner-like, one of kj has different sign than two others. It means that the 
Universe contracts in two directions and expands in the third. The second one, called 
Kasner- unlike, all the parameteres have the same sign, which means that there is a 
contraction in every direction. This is not the end, because from time to time there 
occur short periods in the evolution when spatial derivatives of gravitational field 
dominate over time derivatives, which lead to the transitions between Kasner epochs. 
In each transition the set of the parameters kj characterizing the specific Bianchi I 
model changes to another one. Dynamics of the transitions can be modelled by the 
Bianchi II time evolution [12j . 

In the BKL scheme, the Universe may undergo an infinite number of chaotic-like 
transitions (depending on equation of state of matter fields) from one Kasner epoch 
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to another, and finally collapse to a singularity in a finite proper time [8j [9]. It is clear 
that before the system approaches the singularity, the spacetime curvature acquires 
the Planck scale so the classical scenario cannot be trusted. 

It is clear that, according to the BKL scenario, there are two basic steps in under- 
standing of an early Universe. The first one means the construction of the quantum 
Bianchi I model of the Universe [13], and the second one means finding the quantum 
Bianchi II model. This Thesis is devoted to the first task. 

Quantization in this Thesis is based on the so-called nonstandard Loop Quantum 
Cosmology (LQC) which is an alternative to the standard Loop Quantum Cosmology. 
The difference between them will be discussed later. The standard LQC is a cosmo- 
logical counterpart of Loop Quantum Gravity (LQG) and is is obtained by symmetry 
reduction of LQG, which as a field theory has infinite number of degrees of freedom 
|14j . On the other hand LQG is a non-perturbative theory being today one of the 
most promising candidate to the theory of unification GR with quantum physics. 
Both geometry and matter are dynamical and described by quantum mechanics. In 
LQG there is no background spacetime [13] . 

This Thesis is organized as follows: 
Chapter 1 is an introduction and contains main features of both existing LQC meth- 
ods, namely the standard and the nonstandard. In Section 1 of this chapter we 
describe formalism of the standard LQC concentrating mainly on its successes in re- 
solving the cosmological singularity. Section 2 contains a motivation for another LQC 
method, which underlines the Thesis. In Section 3 we describe this method, called the 
nonstandard LQC and make a comparison between these two methods. In Section 
4 we present an application of the nonstandard LQC to the simplest cosmological 
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model, namely FRW. At the classical level we show the occurrence of Big Bounce. 
At the quantum level particular attention is paid to volume and energy density op- 
erators. Spectra of these operators are analyzed in details. 

Chapter 2 is the main part of the Thesis. Here we present an application of nonstan- 
dard LQC to the Bianchi I cosmological model. Section 1 is devoted to calculations 
done on the classical level, which mean solving the equations of motions and finding 
the algebra of elementary observables. In Section 2 we make, before quantizing, some 
comfortable redefinitions and, what is more important, analyze structure of the phase 
space. So-called true Hamiltonian is introduced. In Section 3 we face the quantization 
in Kasner-like and Kasner-unlike cases. Particular attention is paid to the volume 
operator and the problem of an evolution. 
In Conclusions we make summary of all the results. 

Appendix A presents curvature of connection expressed in terms of holonomies. Clas- 
sification of phase space in terms of observables includes Appendix B. Some ambigu- 
ities in quantization are disscussed in Appendix C. 



Chapter 1 

Loop Quantum Cosmology 



1.1 Formalism of standard LQC 
1.1.1 Basic facts 

By the standard LQC one means LQC that is strongly inspired by LQG [151 HHJ 
[T7] . The inspiration consists mainly in applying the two ingredients of LQG: (i) 
modification of the curvature of connection by loop geometry, and (ii) making use 
of the holonomy-flux variables. The construction of LQC has been carried out by 
mimicry of the construction of LQG. 

The LQC was firstly applied to quantization of FRW model as FRW is the most 
common model in cosmology and astrophysics today. This field is not very old, 
because the first papers appeared less than 10 years ago [TU [19]. Since then more than 
100 papers have been written. Amongst the most important are [HI [201 Ell 1221 123] . 
We recommend also a few review articles [2^1 125[ [26[ [27] . 

The standard LQC [T^l |2"U] means basically the Dirac method of quantization, 
which begins with quantization of the kinematical phase space followed by imposi- 
tion of constraints of the gravitational system in the form of operators acting on a 

5 
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kinematical Hilbert space. Then, one must solve these constraints which means find- 
ing kernels of the operators. The kernels are used to identify the physical Hilbert 
space. 

1.1.2 Big Bounce 

The most important result obtained within the standard LQC is that the classical Big 
Bang (and also the Big Crunch) is replaced by the quantum Big Bounce due to strong 
quantum effects at the Planck scale [2"4"l l2"o] 126] . The result was obtained for FRW 
models with k = and k = 1 and for Bianchi I model [23] . It was done using analytical 
and numerical methods [T£| |2"2"| |2"5] . In all these models the role of internal time was 
played by a scalar field which enables interpreting the Hamiltonian constraint as an 
evolution equation. Singularity is resolved in the sense that observables like energy 
density of matter which classically diverge, are represented by operators bounded 
from above on the states (vectors of the physical Hilbert space) which are semi- 
classical asymptotically [22J . It is suggested in |23j that the bounce may occur for the 
states which are more general than semi-classical at late times, which demonstrates 
robustness of LQC results. Quantum evolution, described by (11.2.240 ; is deterministic 
across the bounce region. An universe undergoes a bounce during the evolution from 
an epoch before the Big Bang to an epoch after the Big Bang, so on the "other side" 
of the bounce there is also a universe [25J. These are main highlights of LQC (see, 
e.g. |27j for a complete list). 

There exists an alternative to the standard LQC approach called the nonstandard 
LQC, which is presented and applied in the rest of the Thesis. 
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1.2 Motivation for nonstandard LQC 
1.2.1 Hamiltonian 

The gravitational part of the classical Hamiltonian, H g , in GR is a linear combination 
of the first-class constraints, and reads [HI EE EE EO [20] 

H g := J d 3 x(N l Q + iV a C a + NC), (1.2.1) 

where £ is the spacelike part of spacetime RxE, (N l ,N a ,N) denote Lagrange 
multipliers, (C», C a , C) are the Gauss, diffeomorphism and scalar constraints. In our 
notation (a, b = 1,2,3) are spatial, and (i,j,k = 1,2,3) internal SU(2) indices. The 
constraints must satisfy a specific algebra. 

For flat FRW model with massless scalar field we can rewrite the gravitational part 
of the classical Hamiltonian, having fixed local gauge and diffeomorphism freedom, 
in the form (see, e.g. [22] ) 

H g = - 7 " 2 / d 3 x Ne~ l e l]k E^E bk F l ab , (1.2.2) 
Jv 

where 7 is the Barbero-Immirzi parameter, V C S is an elementary cell, E is spacelike 
hyper-surface, N denotes the lapse function, e^k is the alternating tensor, E? is a 
densitized vector field, e := y/\ det ~E\, and where F^ b is the curvature of an SU{2) 
connection A l a . 

The resolution of the singularity, obtained within LQC, is based on rewriting 
the curvature F^ b in terms of holonomies around loops. The curvature F% b may be 
determined [22] by making use of the formula (see the Appendix A) 

, h W _ lN 

F k ab = -2 lim Tr ( ° ij ... ) r k °J a (1.2.3) 

aD ArDii-^O V \2T/ 2 / 3 / a a 7 
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where 

h^ = hi%fW)-\h^)- 1 (1-2-4) 
is the holonomy of the gravitational connection around the square loop Dy, considered 

1/3 

over a face of the elementary cell, each of whose sides has length XV with respect 
to the flat fiducial metric °q a b := bij°uj % a °uj^ fiducial triad °e\ and cotriad °oj k a satisfy 
°uj % a °e" = Sp the spatial part of the FRW metric is qw, = a 2 {t)°q ab ; ArD^ denotes 
the area of the square; V a = J v y/°qd 3 x is the fiducial volume of V. Because its value 
is for our analysis not essential, we set Vo = 1. Here we would like to emphasize that 
(I1.2.3P is exact iff A = 0, which is not our case. 

1/3 

The holonomy along straight edge °e a k d a of length \V reads 
r xv 1/3 

h { ^\c) = Texp ( / r (k) A ik) dx a ) = exp(r fe Ac) = cos(Ac/2) I + 2 sin(Ac/2) r k , 
Jo 

(1.2.5) 

where r k = —io^jl [a k are the Pauli spin matrices) and V denotes the path ordering 
symbol. Equation (I1.2.5P presents the holonomy calculated in the fundamental, j = 
1/2, representation of SU{2). 

Making use of ( 11.2.2p . ( 11.2.3P and the so-called Thiemann identity [16] 



s ijk e^E^E bk = |^ °e abc °^Tr{h^{{h^r\V}r) (1.2.6) 

k 



leads to H g in the form 



H g = Jim HV>, (1.2.7) 

A— > 

where 



sgnjp) 
lirGj 3 X 

and where V — \p\% = a 3 is the volume of the elementary cell V. 



H « x) = E w««*i>(aW (htYH^r^uhPr 1 , n), ( i^> 

ijk 
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The connection A\ and the density weighted triad E% which occurs in (I1.2.6P is 
determined by the conjugate variables c and p as follows: A k a = °uj k a c and E% = 
° e k y/QoPi where c = ■ya and \p\ = a 2 . 

It should be noticed that in this section we use the "old" quantization scheme |21j . 
There exists also the "improved" scheme p, = A described in [22J and used in 
section "Application" of this Thesis. It has serious advantages and is now commonly 
used by LQC community. However, obtained results concern both methods [29|. 

The classical total Hamiltonian for FRW universe with a massless scalar field, 0, 
reads 

H = Hg + H^O, (1.2.9) 

where H g is defined by (11. 2. 71) and where sign reminds that if is a constraint of the 
system. The Hamiltonian of the scalar field is known to be: = N p^lp]^^ /2, where 
and p^ are the elementary variables satisfying {(j),p^} = 1. The relation H f» 
defines the physical phase space of considered gravitational system with constraints. 

1.2.2 Quantization 

In the Dirac quantization [301 EI] we find a kernel of the quantum operator H corre- 
sponding to H, i.e. 

HV = 0, (1.2.10) 

since the classical Hamiltonian is a constraint of the system, and try to define a 
scalar product on the space of solutions to (ll.2.10p . This gives a starting point for 
the determination of the physical Hilbert Sp£lC6 1~Lphys' 
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Kinematics 

The classical elementary functions satisfy the relation 

{p,N x } = -i^\N x , (1.2.11) 

where G is the Newton constant. Quantization of the algebra (ll.2.1ip is done by 
making use of the prescription 

{•,•>-£[•,•]• U.2.12) 

The basis of the representation space is chosen to be the set of eigenvectors of the 
momentum operator [13] and is defined by 

p\fi) = —^-fi\fi), (ieR, (1.2.13) 

where I 2 = Gh. The operator corresponding to N x acts as follows 

N x \fi) = |/i + A). (1.2.14) 

The quantum algebra corresponding to (II .2. lip reads 

-[p,N x }\ri = -z—t\NM. (1.2.15) 
in 6 

The carrier space, T g , of the representation (11 .2. 151) is the space spanned by {\fi), fi G 
IR} with the scalar product defined as 

0V):=<W> (1-2.16) 

where S^y denotes the Kronecker delta. 

The completion of J- g in the norm induced by (11.2. 16j) defines the Hilbert space 
T-L g kin = L 2 (M. Bohr ,dfi Bohr ), where M. Bohr is the Bohr compactification of the real line 
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and d\i Bohr denotes the Haar measure on it [H] . l-i 9 kin is the kinematical space of the 
gravitational degrees of freedom. The kinematical Hilbert space of the scalar field is 
Hf. in = L 2 {R, d(p), and the operators corresponding to the elementary variables are 

(#)O) = 0V#), M = -ih^l>. (1.2.17) 

The kinematical Hilbert space of the gravitational field coupled to the scalar field is 
defined to be W, kin = H 9 Un <8> V-tin- 

Dynamics 

The resolution of the singularity [HJ [20j [2TJ [22j [23] is mainly due to the peculiar way 
of defining the quantum operator corresponding to H g . Let us consider this issue in 
more details. 

Using the prescription {•, •} — > •] and specific factor ordering of operators, one 
obtains from d 1 . 2 . 8 [) a quantum operator corresponding to Hg in the form [1] 



6<» -Se«*4T(W)«)-',^)- (12.18) 

P ijk 

One can show [2] that (11.2. 18j) can be rewritten as 

Hi 9 X) M = 8?r7 3 A 3 Z 2 (Va - V) (\» + 4 ^> + 4A» , (1.2.19) 

where is an eigenstate of p defined by ( 11.2.131) . and where is an eigenvalue of 
the volume operator corresponding to V = \p\ 3 ^ 2 which reads 

V\n) = (^^y /2 l 3 p l/i) =: V, \f,). (1.2.20) 

The quantum operator corresponding to H g is defined to be [HJ [21] 

H g := ffW | A=Mo , where < /i„ G 1. (1.2.21) 
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Comparing (ll.2.2ip with (ll.2.7p . and taking into account (I1.2.3P one can see that 
the area of the square is not shrunk to zero, as required in the definition of the 
classical curvature fll.2.3j) . but determined at the finite value of the area. 

The mathematical justification proposed in [IJJ [21] for such regularization is that 
one cannot define the local operator corresponding to the curvature F* b because the 
1-parameter group N\ is not weakly continuous at A = in JF g (dense subspace of 
Hkin)- Thus, the limit A — > of does not exist. To determine \x a one proposes 
in [TJJ |21j [22] the procedure which is equivalent to the following: We find that the 
area of the face of the cell V orthogonal to specific direction is Ar = \p\. Thus the 
eigenvalue problem for the corresponding kinematical operator of an area Ar := \p\, 
due to (11.2. 13p . reads 

Ar = \y\ \n) =: ar(ji) // e K, (1.2.22) 

where ar{fi) denotes the eigenvalue of Ar corresponding to the eigenstate \p). On 
the other hand, it is known that in LQG the kinematical area operator has discrete 
eigenvalues [221 E3] and the smallest nonzero one, called an area gap A, is given by 
A = 2^7^. To identify Ho one postulates in [21] that \i is such that ar(/i D ) = A, 
which leads to \x D = 3v3/2. It is argued [HI EU |22j [23] that one cannot squeeze a 
surface to the zero value due to the existence in the universe of the minimum quantum 
of area. This completes the justification for the choice of the expression defining the 
quantum Hamiltonian (11.2.21I) offered by LQC. 

It is interesting to notice that for the model considered here (defined on one- 
dimensional constant lattice) the existence of the minimum area leads to the reduction 
of the non-separable space T g to its separable subspace. It is so because due to ( 11 .2.1411 
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we have 

iV M » = |// + /i ), (1-2-23) 

which means that the action of this operator does not lead outside of the space 
spanned by + k fi ), k e Z}, where /x G R is fixed. 

Finally, one can show (see, e.g. [HJ [21]) that the equation for quantum dynamics, 
corresponding to (1 1 . 2 . 1 j) . reads 

B(p) ^(/i, 0) - C + (/i)^(/i + 4/i G , 0) - C-{p)^{n - 4/i , 0) - C°(/i)^(/i, 0) = 0, 

(1.2.24) 

where 

W:=(^) 6 0/^ + ^| 3/4 -|^-^| 3/4 ] 6 , C°(ri:=-C + (ri-C-(^ (1.2.25) 

C+(/i) := | |// + 3/i | 3 / 2 - |/i + /i | 3 / 2 | , C-(fi) := C + (fi - 4/i G ). (1.2.26) 

Equation (jl.2.24p has been derived formally by making use of states which belong 
to T := T g ® T<$„ where T g and T§ are dense subspaces of the kinematical Hilbert 
spaces V, 9 kin and T-if, in , respectively. The space T provides an arena for the derivation 
of quantum dynamics. However, the physical states are expected to be in J 7 *, the 
algebraic dual of T (see, e.g. [TJ1 [21] and references therein). It is known that T C 
7-Lkin C . Physical states are expected to have the form < := ^^(/i, 0) < iAi 
where < fi\ is the eigenbras of p. One may give the structure of the Hilbert space 
to some subspace of J 7 * (constructed from solutions to (11.2.24)) ) by making use of 
the group averaging method [MJ ES] and obtain this way the physical Hilbert space 

The argument in 0) is interpreted as an evolution parameter, /j is regarded 
as the physical degree of freedom. Let us examine the role of the parameter \x a 
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in f ll.2.24p . First of all, its presence causes that (ll.2.24p is a difference-differential 
equation so its solution should be examined on a lattice. It is clear that some special 
role must be played by /i Q = as the coefficient functions of the equation, defined 
by (I1.2.25P and (ll.2.26p . are singular there. One can verify [21] that as /i — > the 
equation (ll.2.24p turns into the Wheeler-DeWitt equation 

s d 2 . . IQttG d _9 , . . ArcyGh ,-3/2 

B{ji) Qjpi/)(fi,<p) ——^/jl—'ip(v,(p)=0, with B(p) := | — - — //| 

(1.2.27) 

Equation fll.2.24j) is not specially sensitive to any other value of /x D . Thus, the 
determination of the numerical value of this parameter by making use of the mathe- 
matical structure of (ll.2.24p seems to be impossible. 

1.2.3 Minimum length problem 

The singularity resolution offered by LQC, in the context of flat FRW universe, is a 
striking result. Let us look at the key ingredients of the construction of LQC which 
are responsible for this long awaited result: 

Discussing the mathematical structure of the constraint equation we have found 
that Ho must be a non-zero if we wish to deal with the regular (ll.2.24p instead of the 
singular (11.2.271) . However, the numerical value of fi a cannot be determined from the 
equation (11.2.241) . It plays the role of a free parameter if it is not specified. 

The parameter \x Q enters the formalism due to the representation of the curvature 
of the connection F* b via the holonomy around a loop (11.2.31) . The smaller the loop 
the better approximation we have. The size of the loop, /i , determines the quantum 
operator corresponding to the modified gravitational part of the Hamiltonian (ll.2.2ip . 
One may determine \i by making use of an area of the loop (used in fact as a technical 
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tool). Thus, the spectrum of the quantum operator corresponding to an area operator, 
Ar, seems to be a suitable source of information on the possible values of fi . Previous 
section shows explicitly that the construction of the quantum level is heavily based 
on the kinematical ingredients of the formalism. Thus, it is natural to explore the 
kinematical Ar of LQC However, its spectrum (11.2.221) is continuous so it is useless 
for the determination of On the other hand, the spectrum of kinematical Ar of 
LQG is discrete [321 E3]- Thus, it was tempting to use such a spectrum to fix fi Q 
postulating that the minimum quantum of area defines the minimum area of the loop 
defining (jl.2.2ip . This way /x G has been fixed. 

The physical justification, however, for such procedure is doubtful because LQC 
is not the cosmological sector of LQG. Therefore, Eq. f ll.2.2ip includes an insertion 
by hand of specific properties of the spectrum of Ar from LQG into LQC [35]. After 
all, the area gap of the spectrum of Ar of LQG is not a fundamental constant (like 
the speed of light, Planck's constant, Newton's constant) so its use in the context of 
LQC has poor physical justification. 

1.2.4 Summary 

First of all we have shown that the introduction of the quantum of an area from LQG 
into LQC at kinematical level is only an assumption. As a consequence, the energy 
scale of the Big Bounce described by A parameter is in fact unknown. It is so because 
A is a free parameter of LQC. 
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1.3 Formalism of nonstandard LQC 

Now we will present an alternative method of canonical quantization of cosmological 
models of GR, which makes use of loop geometry (23 EHl EH]- We believe that 
the nonstandard LQC may be related with the reduced phase space quantization of 
Loop Quantum Gravity [ID]. What is the motivation for developing an alternative 
formalism? First of all, agreement of results obtained with both methods would be a 
sort of proof that the procedure of quantization is correct. Of course the final test is 
always an agreement with observational data when they become available. Another 
reason to develop an alternative approach is improving our understanding of some 
conceptual issues like identification of physical observables or quantum evolution of 
a system with the Hamiltonian constraint. 



1.3.1 Main idea 

In the nonstandard LQC [371 138] one first solves the constraint (the constraints) 
at the classical level to identify the physical phase space (i.e. the space of Dirac's 
observables). Secondly, in that space one finds the elementary observables and their 
algebra. These elementary observables are used as "building blocks" for the compound 
observables, like the energy density or the volume of the universe, so they have deep 
physical meaning. The compound observables are thus defined on the physical phase 
space too. Thus, their properties may be confronted in future with the data of 



observational cosmology. T 



re compound observables depend on the elementary ones 



and an evolution parametei 



so for fixed moments of time they are functions only of 



elementary observables. Next step is a quantization. By this term we mean finding a 



which is value of the scalar field 
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self-adjoint representation of the algebra of the elementary observables and solution 
to the eigenvalue problem for operators corresponding to the compound observables 

The difference of understanding the term "quantization" is a source of another dif- 
ference between these two LQC methods. In nonstandard formalism, approximation 
of the curvature of connection by a holonomy along a loop of finite size (modifica- 
tion of Hamiltonian by loops) is done entirely at the classical level. Our approach is 
different from the so-called polymerization method practisized by users of standard 
LQC and treated as an effective quantum theory (see, e.g. [Hj), where the modi- 
fication in the Hamiltonian: /3 — > sin(A/3)/A finishes the procedure of quantization. 
It means that in our method resolution of the singularities happens at the classical 
level due to loop modification of Hamiltonian. This modification is parameterized by 
a continuous parameter A. There is no specific choice of A, so we can say that A is a 
free parameter. Details on theoretical ways of finding the value of A can be found in 
conclusions of this Thesis. 

There is also another important issue. Why should we quantize a cosmological 
model which is free from the cosmological singularity? We have at least three rea- 
sons: (i) to make comparison with the standard LQC results, we must have a quantum 
model; (ii) the parameter A specifying the modification is a free parameter in non- 
standard LQC. As the result, the critical density of matter at the bounce becomes 
unspecified as it depends on A. Since it may become arbitrarily big for small enough 
A, the system may enter an arbitrarily small length scale, where quantum effects can- 
not be ignored [37]; (iii) making predictions of our model for quantum cosmic data 
may be used to fix the free parameter A, after such data become available. 
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1.3.2 Comparison of both LQC methods 

Shortly, one can write: 

the standard LQC = first quantize, then impose constraints = Dirac's method; 

the nonstandard LQC = first solve constraints, then quantize = reduced phase space 

quantization. 

The most important advantage of the nonstandard LQC is that the spectra of the 
operators are directly obtained on physical Hilbert space. In the case of the standard 
LQC, one firstly obtains results on the kinematical Hilbert space. The physical states 
are obtained using kernels of quantum constraints operators. Applying group averag- 
ing metods leads to physical spectra of observables. Thus our nonstandard method 
is simpler and more efficient than the standard one. 

Another important feature of nonstandard LQC is that this method is fully con- 
trolled analytically as it does not require any numerical work, at least in FRW and 
Bianchi I cases, in contrast to the standard LQC results. 

In the nonstandard LQC an evolution parameter stays classical during the quan- 
tization. This happens because does not belong to the physical phase space. In the 
standard LQC is a phase space variable and should be quantized. This is crucial be- 
cause being a quantum variable may fluctuate, which may makes an interpretation 
of problematic. 

1.4 Application of nonstandard LQC to FRW 

In this chapter we consider, as an application of our method, the simplest cosmological 
model, namely FRW model with k = and with free massless scalar field in space 
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with topology 



1.4.1 Hamiltonian 

In what follows we use the "improved" scheme p, = ■sj^ A [22]. Using it and taking 
( I1.2.5P we calculate (I1.2.8P and get the modified total Hamiltonian corresponding 
to (THT9]) in the form 

H M/N= ^ Sin2(A/3) , + ^ (141) 

1 8vrG 7 2 A 2 + 2v' 1 J 

where 

P'-=^> -=N 3/2 (1-4.2) 

are the canonical variables proposed in [22] ■ The variable = ja/a so it corresponds 
to the Hubble parameter a/a, whereas t> 1//3 = a is proportional to the scale factor a. 
The complete Poisson bracket for the canonical variables (/3, v, <p,p<f,) is defined to 

be 

d- d- d- d- 



{•,-}:=4ttG 7 



d- d- d- d- ^ 



d(3 dv dv d(3 
The dynamics of a canonical variable £ is defined by 

£:={£,H^}, ^e{/3,v,<f>,p^}, (1.4.4) 

where £ := d^/dr, and where r is an evolution parameter. The dynamics in the 
physical phase space, J^ ys , is defined by solutions to (11 .4. 41) satisfying the condition 
ffW « 0. The solutions of (jTOl ignoring the constraint if( A ) « are in the 
kinematical phase space, J 7 ^. 
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1.4.2 Classical dynamics 

Equation (11.4. ip can be rewritten as 

H w = NH [ X) H {X) « 0, (1.4.5) 

where 

where k 2 = 47rG/3. 

It is clear that Hq^ = only in the case when = = sin(A/3). Such case, due to 
(ll.4.7p - (ll.4.1ip . implies no dynamics. 

Choosing the gauge N := (which simplifies the calculations) we get 

P4> = 0, (1.4.7) 

$ = -4vrG 7 ^M, (1.4.8) 
A 

= K7 sgn(^), (1.4.9) 

v = AnG-fv cos(\ (3) sgn(sin(A /3)), (1.4.10) 

H w = 0. (1.4.11) 

Combining fll.4.9j) with (11.4. 10p gives 

v 

— = 3kv cos (A/3) sgn(sin (A/3)) sgn^). (1.4.12) 



Rewriting (11.4.121) (and using v/<p = dv/dcj)) gives 

sgn(sin(A/3)) dv , u , . 1Q x 

__ = 3Ksgn(p^)# 1.4.13 

cos(Ap) v 

Making use of the identity sin 2 (A/3) + cos 2 (A/3) = 1 and (11.4.111) gives 



ms(m = ^.^y (1 .4.i4) 
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Combining ( 11.4.130 with fll.4.14j) . for /3 g]0,tt/2A[, leads to 

dv 



3k sgnQo^) c£0. (1.4.15) 



3«sgn(p^)(0-0 o ). (1-4-16) 



V^ 2 - («7^P 

Since is just a constant (due to (11.4.70 ) we can easily integrate (11.4.150 and get 

In v + \Jv 2 - (i^l^P<f>y 
Rewriting (11.4.161) leads to 

2v = exp (3k sgn(^) (0 - O )) + («7b*l • exp ( - 3k sgn(p ) (0 - O ))- (1-4.17) 
The solution for the variable [3 may be easily determined from (11.4.110 rewritten as 

K7W = „^M (1 .4.18) 

Finally we get 

2k7A|^| 



sin(A/3) 



exp (3«sgn(zty) (0 - O )) + (K7MO 2 exp ( - 3/csgn(p ) (0 - <f> )) 

(1.4.19) 



where the domain of the variable /3 has been extended to the interval ]0, 7r/A[. 

Equations (11.4.170 and (11.4.190 present the dependence of the canonical variables 
v and (3 on the evolution parameter 0, which is a monotonic function due to (11.4.90 . 

1.4.3 Observables 

Elementary observables and their algebra 

A function, O, defined on phase space is a Dirac observable if 

{0,H {X) } « 0. (1.4.20) 
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Since we have 

{0,H W } = {O,NH ( X] H w } = NH ( X) {O,H w } + {£>, NH ( X) }H {X) , (1.4.21) 

it is clear that on the constraint surface, = 0, the Dirac observable satisfies 

(independently on the choice of N) a much simpler equation 

{O,H (x) }^0. (1.4.22) 

Using the gauge N := l/H^ and solve f ll.4.20p in the whole phase space, i.e. we 
solve the equation 

sm(A/3) dO DO Ksgn(p ) dO 

— - v cos(A/3) — — — = 0. (1.4.23) 

A function O = 0{O u ...O k ) satisfies fll.4.231) if 

{0 1: ffW} = = {0 2 , H^} = ... = {O k , H^}, (1.4.24) 

where k + 1 is the dimension of the kinematical phase space. It is so because one has 

{O, H™} = ^ {0 1} ffW} + . . . + £L {O kl H^}. (1.4.25) 

In what follows we consider only elementary observables. The set of such observ- 
ables, £, is defined by the requirements: (i) each element of £ is a solution to f l 1 . 4 . 2 3 j) . 

(ii) elements of £ are functionally independent on the constraint surface, H^ x ' = 0, 

(iii) elements of £ satisfy a Lie algebra, and (iv) two sets of observables satisfying two 
algebras are considered to be the same if these algebras are isomorphic. 

In our case k = 3 and solutions to fll.4.231) are found to be 

0i:=j>0, 2 :=ct>-^- arth(cos(A/3)), 3 :=sv^^-, (1.4.26) 
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where s := sgn(p^). One may verify that the observables satisfy the Lie algebra 

{0 2 ,0 1 } = 1, {O x ,O 3 } = 0, {0 2 ,0 3 } = y K . (1.4.27) 

Because of the constraint = (see (I1.4.18P ). we have 

C» 3 = 7 kC 1 . (1.4.28) 

Thus, we have only two elementary Dirac observables which may be used to param- 
eterize the physical phase space Fph ys - To identify the Poisson bracket in J^ ya con- 
sistent with the Poisson bracket (I1.4.3P defined in J 7 ^, we find a symplectic twoform 
corresponding to (11.4. 3p . It reads 

u = 1 d/3 A dv + dtp A dp*. (1.4.29) 
47rG7 

The twoform u is degenerate on J 7 ^,, due to the constraint H ^ = 0. Making use 
of the explicit form of this constraint ( 11.4.181) and the functional form of 0\ and 0%, 
leads to the symplectic form f2 on J 7 ^,,,. Direct calculations give 

: = U\hw=o = d0 2 AdO u (1.4.30) 

where u<g(\) =0 denotes the reduction of uj to the constraint surface. The Poisson 
bracket corresponding to (I1.4.30P reads 

, , d- d- d- d- 

:= dO~ 2 dO~i~ dO~idO~2 (1-4.31) 

so the algebra satisfied by 0\ and O2 has a simple form given by 

{0 2 ,0 1 } = 1. (1.4.32) 

Our kinematical phase space, J^j^l, is four dimensional. In relative dynamics 
one variable is used to parameterize three others. Since the constraint relates the 
variables, we have only two independent variables. This is the reason we have only 
two elementary physical observables parameterizing J 7 ^^. 
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Functions on phase space 

Now we discuss the functions on the constraint surface that may describe singularity 
aspects of our cosmological model. Considered functions are functions of elementary 
observables and evolution parameter 0, so they are not observables. They become 
observables for fixed cf). 

An interesting function is energy density p of the scalar field <fi 

,2 



p(A,0) = ^. (1.4.33) 

— u 



In terms of elementary observables we have 

p <j> = Oi, v = n-f\\0 1 \ cosh (3k(0 - £> 2 )) (1.4.34) 

which means that 

p(X,<f>) = - 1 . (1.4.35) 

HK J 2 (n-fX) 2 cosh 2 3n((j) -0 2 ) 

For fixed p^ the density p takes its maximum value at the minimum value of v. 
Rewriting (11 .4. 171) in the form 

v 

— = cosh (3ks(0 — 4>o) — In A), where A := firyA (1.4.36) 

we can see that cosh(-) takes minimum value equal to one at 3ns (<j) — <p ) = In A. 
Thus, the maximum value of the density, p max; corresponds to v — A and reads 

Pmax = Tjj ■ (1.4.37) 

which means that the Big Bounce occurs at the classical level due to modification 
classical Hamiltonian by loops. We can determine p max if we know A, but as we now 
A is a free parameter of the formalism. 
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Another interesting observable is a volume. In terms of observables it reads [37] 
v(X,(f>) = KjXpx] cosh3/t(0- 2 ). (1.4.38) 

Looking at the equations written above one can see that p and V depend explicitly 
on observables 0\ and 2 and evolution parameter so they are observables for fixed 
value of 0. 

1.4.4 Quantization 

As we remember, apart from resolving the singularity there are some important rea- 
sons to continue the procedure of the quantization. 

Representation of elementary observables 

For the classical algebra ( jl.4.32[) two representations are used 

1 — > 3 1 f(x) := -ihd x f(x), 2 — ► 6 2 f{x) :=xf(x) :=xf(x), (1.4.39) 
which leads to [Oi, 2 ] = —ihl, and 

1 — ► dj(x) := xf(x) := xf(x), 2 — > d 2 f(x) := -ihd x f(x), (1.4.40) 
which leads to [Oi, 2 \ — i hi, where iGl 

Due to the Stone— von Neumann theorem all self-adjoint representations of the algebra 
(11.4.32p are unitarily equivalent to the representation ( 11 .4.391) or (11.4.40j) defined on a 
suitable dense subspace of L 2 (M). In that sense the choice of representation is unique. 
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Energy density operator 

The representation (ll.4.40p is essentially self-adjoint on the dense subspace D of the 
Hilbert space L 2 [—r,r], where r e R + , defined to be 

D:={f e C°°[-r, r] | f {n) {-r) = f {n) {r), n G {0} U N}, (1.4.41) 

where := d n f/dx n . 

The eigenvalue problem, 02 /p = p/ P , has the solution 

/p(x) = (2r)~ 1/2 exp(ia;p//i), p(Jfe) := 2vr^/r, fceZ. (1.4.42) 

The spacing of neighboring eigenvalues □ is defined by 

□ :=p(k + 1) -p(k) = 2^h/r (1.4.43) 

Making r sufficiently large □ can be made as small as desired, which means that the 
spectrum of O2 is continuous. 

In the representation (jl.4.40p the energy density operator reads 

p:=- a - • (1.4.44) 

2(/t7A) 2 cosh 3k(0 + i hd x ) 

Since 2 is essentially self-adjoint on T r := {f p (k)}kez, we may apply the spectral 
theorem to get pf p = p(<p, X,p) f p where 

^ ^' 2 (K7A) 2 cosh 2 3k(0 — p) ' ^ 

and where p(0, A,p) is the eigenvalue corresponding to the eigenvector f p . 

It is clear from our results that classical (11.4.331) and quantum (11.4.451) expressions 

for the energy density coincide. One may verify that the maximum density p m ax{^) = 
1 1 



2 (k 7 A) 2 • 

Starting from the other representation we would get the quantum model of the 
energy density presented in [J2], which is equivalent this present one. 
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Volume operator 

To define the quantum operator corresponding to v, we introduce w defined by 

w := KjX O x cosh3/c(0 - 2 ). (1.4.46) 

Since v — \w\ it is clear that quantization of v reduces to the quantization of w. The 
latter may be done in a standard way as follows 

wf(x) := kjX ~ (di cosh 3/c(0 - <5 2 ) + cosh 3/c(0 - 2 ) (1.4.47) 

where / G L 2 (M). 

For the elementary observables 0i and O2 we use the representation (11.4.39[) . An 
explicit form of the operator w reads 

w f(x) = i — ^ — ^2 cosh 3k(0 — x) 3k sinh 3/t(0 — a;)^ f(x). (1.4.48) 

To simplify our considerations we take / in the form 

f(x) := A e ih{x) cosh~ 1/2 3re(</> - x), (1.4.49) 

where h is a real- valued function and AeK. 

Eigenvalue problem 

Considering the eigenvalue problem for the operator w in the set of functions of the 
form fll.4.491) . We get 

wf(x) = -^XHcosh3K((f) - x)^^- f(x) =: bf(x), (1.4.50) 

ax 

where b G M. is the eigenvalue of w. 
A general form of h satisfying f l 1 . 4 . 5 j) is given by 

2b 

h(x) = - — arctan e ^- x \ (1.4.51) 
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which means that a normalized /& satisfying (jl.4.50p reads 

3 3k(</>— x) 



h{x) ■-- 



' K exp ( i 3 J^ X n arctan < ' 



cosh 2 3k( 



x 



Orthogonality 



Using (I1.4.52P we get 



(Mfa) = - 

7T 



3k r ex Pl z S3 arctane3K( "" x) 



, cosh3fi;(0 — x) 

The substitution tan z = e 3K ^~ x ' leads to 



dx. 



One may verify that {fb\fa) = iff 



3k 2 7A/i / 2(a-b) 
-% — — exp I % -r-^. — —z 



ir(a — b) 



a — b = 6K 2r y\hm = SnG'jXhm, mGZ. 



;i.4.52) 



(1.4.53) 



1.4.54) 



1.4.55) 



Thus, the set T b :— { f a \ a — b + 8irG'-f\hm] m G Z; b G R } is orthonor- 
mal. Each subspace Tb C L 2 (R) spans a pre-Hilbert space. The completion of each 
span J 7 !,, Vfe G R, gives L 2 (R) in the norm of L 2 (R). 



Self-adjointness 

The operator w is symmetric on T b for any b G R due to (fb\wf a ) — (wfb\f a ) = 
(a - b)(f b \fa) because (f b \f a ) = for a ^ b. 

To examine the self-adjointness of the unbounded operator w, we first identify the 
deficiency subspaces of this operator K,± |I3J 



£± := {g± G D b (w*) | (^±|(u>±il)/ a ) = 0, V/ G A,(t/>)}, 



1.4.56) 
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where D b (w) := span JF h , and D b (w*) := {/ G L 2 (R) : 3!/* = (f\wg), Vg G 

For each f a G D b (w) C L 2 (M) we have 

/oo 
dx~g~dx)f a (x) =► <? + = = <7_. (1.4.57) 
-oo 

Thus, the deficiency indices n± := ofam[/C±] of u> satisfy the relation: n + = = 
This proves essential self-adjointness of u> on D b (w). 

Spectrum 

Due to the spectral theorem on self-adjoint operators [43l|44], we may quantize the 
volume as follows 

v =\ w \ — ->. f)/ a : = | a |/ a . (1.4.58) 

A common feature of all is the existence of the minimum gap A := STcG'jhX 
defining a quantum of the volume. Let us discuss this issue in more detail. Denoting 
the minimum eigenvalue of v by v m i n , one can verify that v m i n = min{6, A — b}, where 
b G [0, A[. The spectrum consists of the union of {v m i n + nA} and {— v m i n +(n+l)A}, 
where n = 0, 1, . . . There are only two cases when these two subsets are identical, 
namely when v m %n = or v m { n = A/2 , for which the minimum gap A is a constant 
gap between any two adjacent levels of the spectrum. Otherwise, the gap equals 
either A — 2v min or 2t> min , and the minimum gap is the smaller one. One can verify 
that the case of any b G K. reduces to the above case. 

There is no quantum of the volume in the limit A — > 0, corresponding to the 
classical FRW model without the loop geometry modification. 

It results from ( I1.4.55P that for 6 = and m = the minimum eigenvalue of 
v equals zero. This special case corresponds to the classical situation when v = 0, 
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which due to (jl.4.ip means that = (no classical dynamics). Thus, we have 
a direct correspondence between classical and quantum levels corresponding to this 
very special state. All other states describe bouncing dynamics. 

Evolution 

It is clear that the relation between eigenvectors corresponding to the same eigenvalue 
for different values of the parameter reads ft + ^ = e^ 9 ^ ft — e ~ 1 ^° 1 ft ■ 
One may verify that 

M + 4j) = cosh (3«^W0) + smh ( 3/t ^) g tibU), (1.4.59) 

6K 

thus 

VtW + VO/f) = {ft\<M)fi) cosh (3«0) + smh ^\ fi\dM<P)f^) 

Orb 

= a cosh (3i#) S ab + (b- a f inh ^\ ft\d x fi). (1.4.60) 

OK 

An evolution of the expectation value of the operator w is found to be 

(f(((>)\w((f> + ^)f ((/>)) = Acos\i3k(iIj + B), (1.4.61) 

where / := J2 a afa, fa G ^6- 
One may verify that 

1 , X + Y 



A = sgn(X)VX*-YZ, B = —\n —— , (1.4.62) 

ok X — Y 



where 



v-^i .2 t r "fe«a - OL a a h m(2a + <omK 2 ^h\) . no\ 

X := > a a a, r := > ; — -, (1.4.63) 

^' 1 ' ^ in (2m- l)(2m + l) ' V ; 

a a, m \ / \ / 

and where b = a + 6K 2, yXh, 6 e 1, m G Z, and |X| > |y|. 

One can see that the evolution of the expectation value of the operator w coincides 
with the classical expression fll.4.46[) . 
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1.4.5 Summary 

The resolution of the cosmological singularity is due to the loop modification of the 
Hamiltonian already at the classical level. This modification is parameterized by a 
continuous parameter A, which value is so far unknown. Each value of that parame- 
ter specifies quantum of the volume and the maximum energy density of the matter 
(scalar field in our model). 

Spectrum of the the volume operator is bounded from below and discrete. Its expec- 
tation value coincides with the classical expression. 

The spectrum of the energy density is bounded from below and continuous. There 
is a coincidence between classical and quantum expressions for the maximum energy 
density. 



Chapter 2 



Bianchi I model in terms of 
nonstandard LQC 



From Introduction we know that the Bianchi I model of the Universe is of primary 
importance as it underlies, to some extent, the Belinskii-Khalatnikov-Lifshitz (BKL) 
scenario [TTj H5J HEl H7J HE], which is believed to describe the Universe in the vicinity 
of the cosmological singularity It has been examined recently within the nonstandard 
LQC PS, [50], and has been also studied in the context of the standard LQC [T3J, [191 
[5T1 152| 153] . Clear exposition of the singularity aspects of the Bianchi I model can be 
found in [SI [48]. 

This chapter presents analyzes of the Bianchi I model within the nonstandard 
LQC formalism. In what follows we consider the Bianchi I model with free massless 
scalar field in T 3 topology. We choose this topology, because taking M 3 might lead to 
problems with interpretation of the spectra of the volume operator [50J. 
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2.1 Classical level 
2.1.1 Hamiltonian 

The Bianchi I model with massless scalar field is described by the metric: 

3 

ds 2 = -N 2 dt 2 + J2 a2 i(t) dx l (2.1.1) 

i=l 

where 

/ t \ kl 3 3 

fli(r) = fli(O) ( - J , dr = Ndt, ^^ = 1 = Jjk' + kJ, (2.1.2) 

V r °/ i=l i=l 

and where describes matter field density (k^, = corresponds to the Kasner model). 

To make this chapter self-contained we remind that in GR the gravitational part 
of the classical Hamiltonian, H g , is a linear combination of the first-class constraints 

H g := [ d*x{N i C i + N a C a + NC), (2.1.3) 

where S is the spacelike part of spacetime KxE, (N l , N a , N) denote Lagrange 
multipliers, (Cj, C a , C) are the Gauss, diffeomorphism and scalar constraint functions. 
In our notation (a, b = 1,2,3) are spatial and (i,j,k = 1,2,3) are internal SU(2) 
indices. As we know the constraints must satisfy a specific algebra. 

Having fixed local gauge and diffeomorphism freedom we can rewrite it for the 
Bianchi I model with massless scalar field jl9] 

H g = -r 2 [ d 3 x Ne- 1 e ijk E a °E bk F* b , (2.1.4) 
Jv 

where 7 is the Barbero-Immirzi parameter, V C X is an elementary cell, S is spacelike 
hypersurface, iV denotes the lapse function, is the alternating tensor, Ef is a 
densitized vector field, e := yj det E\, and where F* b is the curvature of an SU(2) 
connection A l a . 



34 



As we know, the resolution of the singularity, obtained within LQC, is based on 
rewriting the curvature F* b in terms of holonomies around loops. The curvature F* b 
may be determined by making use of the formula 

where 

h n .. = h^h^ih^)- 1 ^)- 1 (2.1.6) 

is the holonomy of the gravitational connection around the square loop Dy, considered 
over a face of the elementary cell, each of whose sides has length fijLj (and V a := 
LiL 2 L 3 ) with respect to the flat fiducial metric °q a b '■= ^ij 0(jJ a 0(jJ ai the fiducial triad 
°e£ and cotriad °cu^ satisfy °oj % a °e a - = Sy, Ar Oy denotes the area of the square; and 
= Jy y/°qd 3 x is the fiducial volume of V. 

The holonomy in the fundamental, j = 1/2, representation of £77(2) reads 

h ( t ] = cosGu^/2) I + 2 sinOiQ/2) n , (2.1.7) 

where T{ = —iai/2 (crj are the Pauli spin matrices). The connection A\ and the 
density weighted triad Ef (which occurs in ( 12.1. lip ) are determined by the conjugate 
variables c and p: 

A\ = Cl L- 1 X, E« = Pl Lf 1 L k -' °el (2.1.8) 



where: 



and 



c i = ^d i L il \Pi\ = djduLj L k . (2.1.9) 



{a,p j } = SKGhf6 ij (2.1.10) 
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Making use of (I2.1.4p . (" 12 . 1 . 5 [) and the so-called Thiemann identity 

-1 vai vbk Sgn(piP2p 3 ) 



S E °^Tr^^TW}r) (2,,!) 



leads to H g in the form 



# 9 = lim n (2.1.12) 

where 



Mi>M2,M3-s > 



^»») = - s g"^ lP2P3) Y Ne^Trfh^h^ih^)- 1 ^)- 1 

y 21T&Y /ii/i 2 /U 3 ^-f V J J 

(2.1.13) 



and where V = a\ a 2 a 3 LiL 2 L 3 is the volume of the elementary cell V. 

The total Hamiltonian for Bianchi I universe with a massless scalar field, 0, reads 



H = H g + H <p ^0, (2.1.14) 

where H g is defined by (12. 1. 12[) . The Hamiltonian of the scalar field is known to 
be: = N pV^p^p^ 2/2, where and are the elementary variables satisfying 
{4>,P<f>} = 1- The relation H ps defines the physical phase space of considered 
gravitational system with constraints. 

Making use of (12 . 1 . T[) we calculate ( 12.1.13)) and get the modified total Hamiltonian 
Hg X ^ corresponding to (12.1. 14[) in the form 



H [x) /N = 1 sgn(PiP2P 3 ) 



sin(ci/ii) sin(c 2 /z 2 ) ^3 sgn(p 3 )W PlP f +cyclic 



|P3 



+ 



2V 
(2.1.15) 



where 



J^jA, (2.1.16) 
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and where A is a regularization parameter. Here we wish to emphasize that (12.1.151) 
presents a loop modified but classical Hamiltonian. 

It is known [13], HU [52] that such a choice of leads to the dependence of the final 
results on the fiducial volume Vq. In the universe with compact topology, like that 
considered by us, Vq has physical sense, in contrast to a case with noncompact topol- 
ogy. Thus, an expected dependence of the results on Vq would be rather meritorious 
than problematic. 



In the gauge N = P2 P3~| the Hamiltonian modified by loop geometry reads 
1 



\PiP2\ 3/2 sin(ci/ii) sin(c 2/ U2) + cyclic 



8ttG 7 2 A 2 

The Poisson bracket is defined to be 

3 



(2.1.17) 



{-,-}:= 87tG 7 



k=i 



d- d- d- d- 



dc k dp k dp k dc k 



+ 



d- d- d- d- 



(2.1.18) 



where (ci, C2, cs,pi,p2,P3, 4>,P<f>) are canonical variables. The dynamics of £ reads 



i ■= {^H (X) }, € e {ci,C2,C 3 ,Pl,P2,P3,((>,P<t>}- 



(2.1.19) 



The dynamics in the physical phase space, J^L, is defined by solutions to ( 12.1. 19ft 



satisfying the condition rj 0. The solutions of (12. 1 . 19j) ignoring the constraint 
Rj are in the kinematical phase space, J 7 ^. 



We use the following canonical variables 



\Pi 



Vi ■= \Pi\ 3/2 , 



(2.1.20) 



where % = 1,2, 3. They satisfy the algebra 



{favj} = 12irGy5. 



(2.1.21) 
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where the Poisson bracket reads 

3 



{•,-} = 127rG 7 J2 



k=l 



d- d- d- d- 



<9ft dv k dv k <9ft 



d- d- d-d- , _ 
. (2.1.22 



The Hamiltonian in the variables (12.1. 2U|) turns out to be 



V% 1 /sin(Aft) sin(Aft) sin ( Aft ) sin (Aft) 



ti = — - _ 9 ri5 u l v 2 + 



2 87rG7 2 V A 2 A 

sin(Aft) sin(Aft) 
"A 



f «uu^ a ,p^ 8 ; I (2123) 



where A parameterizes the holonomy of connection modifying the Bianchi I model. 

2.1.2 Classical dynamics 

The Hamilton equations of motion read 

ft = -187rG^^(0 J + O fc ), (2.1.24) 

Vi = 18vrG^cos(Aft) (O,- + O fc ), (2.1.25) 

= P<j>, (2.1.26) 

= 0, (2.1.27) 

H w « 0, (2.1.28) 

where z, j, = 1, 2, 3 and i ^ j ^ k and where 

„ Vi sin (Aft) 

Oi := ^-V^. 2.1.29 

12vrG7A V ; 

Insertion of f!2.1.25p into 02.1.24p gives 

dft = -^M^, (2.1.30) 
A Vi 



which leads to 

sin 

1): 

A 



sin(Aft) 

fi = const (2.1.31) 
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Therefore, Oj are constants of motion. 



Making use of f!2.1.26p . 02.1.25P and cos(A/3;) = ^1 - sin(AA) 2 gives 



dvi 



18ttG 



(Q 3 - + Qa 

P<t> 



Integration of (I2.1.32p leads to 
In 

Thus we have 



Vi + Jv? - (12ttG 7 A Oi) 



18ttG 

P<j> 



(0, + A 



18ttG 



2 \vi\ = exp (O,- + O k ) (0 -</>•) + (127rG7AOi) 2 x (2.1.34) 



x (> X p ( - 1^— (Oj + O fc ) (0 - 0°) , . 



which may be rewritten as 

v< = 127tG 7 A |0<| cosh (1^5 (Qj. + O fe ) 



(2.1.32) 



(2.1.33) 



^)-ln|l27rG?7AO<|J. (2.1.35) 



2.1.3 Elementary observables 

As we know, function F defined on the phase space is a Dirac observable if it is a 
solution to the equation 

{F,H w }^0. (2.1.36) 
An explicit form of (I2.1.36[) is given by 



12ttG 7 ^( 

i=i \ 



; 'dFdHW 8FdHW\ dF 

+ — 



d/3i dvi dvi d/3i J 



which reads 



18ttG 



t>i cos (A/% 



dF_ _ sin(Aft) dF 
dvi A 



<9F 



(2.1.37) 



(O,- + O fc ) + — ^ = 0. (2.1.3* 
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Kinematical observables 



One may easily verify that Oj satisfy (12.1.38ft . Instead of solving ( I2.1.38[) one may 
use the constants that occur in (12.1.351) . This way we get 



Ai = In 



tan (Ml) 



18ttG 



(O, + A 



(2.1.39) 



The observables (12. 1.39ft are called kinematical as they are not required to satisfy 
the constraint fl2.1.28|) . 



Dynamical observables 



An explicit form of the constraint (12. 1.28ft in terms of Oj is given by 



sgn(^) = 6v^G VOxOs + 0^0 3 + 2 0, 



(2.1.40) 



It results from (22), (27) and (28) that dOa + Oi0 3 + 2 3 > so (42) is well 
defined. Thus, the dynamical observables, A^ yn , corresponding to (12.1.39ft read 



A? yn = In 



tan(M) 



+ 



3VttG sgn(^) (Oj + O k ) 
VO1O2 + O1O3 + O2O3 



Algebra of elementary observables 

One may verify that Af yn satisfy the following Lie algebra 

{0,,0,} = 0, 
{Af n ,0,} = 

{A dyn )A dyn } = Q 



(2.1.41) 



(2.1.42) 
(2.1.43) 
(2.1.44) 



In the physical phase space the Poisson bracket is found to be 



{-,} 



dyn 



i=l 



d- d- 



d- d- 



(2.1.45) 
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and the algebra reads 

{0^,0, } dyn = 0, (2.1.46) 

{Af n ,0,} dyn = dy, (2.1.47) 

{Af n ,Af n } dyn = 0. (2.1.48) 

2.1.4 Compound observables 

In what follows we consider the physical observables which characterize the singularity 
aspects of the Bianchi I model. It is helpful to rewrite f)2.1.40p and (12.1.351) in the 
form 

pi = 36ttG (Oi0 2 + dOs + 2 3 ), (2.1.49) 

v t = 12vrG 7 A | 0,| cosh ( ^ & + ° k) + In ( *) - A**") . (2.1.50) 

1 11 V VO!0 2 + 0x03 + 0203 \2) 1 ) 1 ' 



Pi(A»:=^. (2.1.51) 



The so-called directional energy density |4T] is defined to be 

2v\ 

The bounce in the z-th direction occurs when pi approaches its maximum 
which happens at the minimum of Vi (p^ is a constant of motion). One may easily 
verify that in the case when all three directions coincide, which corresponds to the 
FRW model, these densities turn into the energy density of the flat FRW with massless 
scalar field [37]. 

It is clear that takes minimum for cosh(-) = 1 so we have 

„r» = 12 lG7 AO„ ^ = I( TS ^-) 2 . (2.1.52) 



Rewriting Oj and p^ in terms of kj and (41 



O, = | k, K, p = V^Gk^K, (2.1.53) 
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where K is a constant, leads to 



167rG 7 2 A 2 \ki 



) 



2 



ft 



max 



(2.1.54) 



We can determine p™ ax if we know A, but as we remember A is a free parameter of 
the formalism. 

One may apply (I2.1.54P to the Planck scale. Substituting X = Ipi gives 



which demonstrates that p^ ax may fit the Planck scale depending on the ratio k^/kj. 



Another important physical observable is the volume of the Universe. From the 
definitions ( jgXgD and f!2.1.20p we get 



It is clear from ( 12. 1.501) . (12. 1.53ft and (12. 1.21) that the volume is bounded from below. 

2.2 Preparations to quantization 

This section is devoted to some redefinitions which are helpfull to the procedure of 
a quantization. The quantization is required despite the fact that the singularity 
problem is resolved already at the classical level due to the modifications based on 
the loop geometry. The reasons are described above. Here we analyze the structure 
of the phase space, which is obviously much more complicated than in the FRW case. 




(2.1.55) 



V = a x a 2 a z = (v 1 v 2 v 3 ) 



1/3 



(2.1.56) 
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2.2.1 Redefinitions 

Redefinition of evolution parameter 

Firstly we slightly redefine the elementary Bianchi observables used in previous sec- 
tions. Now they read 

_ ;= i ujjtajAft) 

3^7 A 

and 

A - 1 lnf l taJ1 ^h I 3 s 8*M(0j + 0k)<i> (2 2 2) 

" 3k V \ / 2 v / 3v / OlO^T0^07TO^O;' 

where n 2 := 4tiG/3. One may verify that the algebra of redefined observables is 
isomorphic to the previous one 

{O i ,O,-} = 0, {^,0,}=^-, {A l ,A,} = 0. (2.2.3) 

and 

, , / 3v^G (0 7 + O fc ) . \ , 

Since the observables Oj are constants of motion in e R, it is possible to make 

the following redefinition of an evolution parameter 

v3 6 _ n . 
V : = — y 2.2.5 
2 v / 1 2 + Oi0 2 + 2 3 

so we have 

Ui = 3«7A|Oj| cosh3/t((O i + O fc ) <p - A;), (2.2.6) 
which simplifies further considerations. 

New elementary observables 

One can make the following redefinitions 

Ai := Ai - (Oj + O fc ) (p. (2.2.7) 
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Thus, the directional volume (I2.2.6P becomes 



Vi := \Wi 



Wi = 3/t7AOj cosh(3nAi). 



(2.2.8) 



The algebra of observables reads 



{O,,O,} = 0, {A l ,0 J } = 5 t 



\^4.j, Aj } 0, 



(2.2.9) 



where the Poisson bracket is defined to be 




dO k dA k 



d- d- 



) 



(2.2.10) 



2.2.2 Structure of phase space 

All considerations carried out in the previous section have been done under the as- 
sumption that the observables Oi, O2 and O3 have no restrictions. The inspection 
of ( 12.2.21) . ( 12.2.41) and ( 12.2.71) shows that the domain of definition of the elementary 
observables reads 



where k = 1,2, 3. The restriction O1O2 + O1O3 + O2O3 > is a consequence of the 
Hamiltonian constraint (see, [49] for more details). 
In what follows we consider two cases: 

1. Kasner-unlike dynamics: (a) Oj > 0, Oj > 0, Ok > 0, which describes all three 
directions expanding (b) < 0, 0^ < 0, < 0, with all directions shrinking. 

2. Kasner-like dynamics: (a) Oj > 0, > 0, < 0, which describes two 
directions expanding and one direction shrinking; (b) Oj < 0, < 0, Ofc > 0, 
with two directions shrinking and one expanding. 



D := {(A k , O k ) \A k eR, 0,02 + 0,0s + 2 3 > 0}, 



(2.2.11) 
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This classification presents all possible nontrivial cases. Our terminology fits the one 
used in [51] due to the relation Oj = 6K,kiK, (0 < K = const), where constants 
are defined by f[2~l~2]) . 

For more details see Appendix B. 



2.2.3 True Hamiltonian 

Now we can define a generator of an evolution called a true Hamiltonian EL Making 
use of f )2.2.7p . and Oj = const (see [49J), we get 

{A, M} := ^ = -(O, + O fc ), {Oi, H} := ^ = 0. (2.2.12) 

d(p dip 

The solution to (12.2.121) is easily found to be 

EI = Oi0 2 + Oi0 3 + 2 3 . (2.2.13) 

The true Hamiltonian is defined on the reduced phase space which is devoid of con- 
straints. It generates a flow in the family of volume quantities, enumerated by the 
evolution parameter. 



2.3 Quantum level 

2.3.1 Representation of elementary observables 

We use the Schrodinger representation for the algebra f!2.2.9|) defined as 

- h d - 

O k ^ O k fk(x k ) := - — f k (x k ), A k -)> Akfk(x k ) := Xkfk(xk), k = 1,2,3. 

(2.3.1) 
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One may verify that 

[Oi, Oj] = 0, [X, Aj] = 0, [X, Oj] = ifrSij- (2.3.2) 

The representation is defined formally on some dense subspaces of a Hilbert space to 
be specified later. 

2.3.2 Kasner-unlike case 

The condition O1O2 + O1O3 + O2O3 > is automatically satisfied in this case, because 
Oi, O2 and O3 are of the same sign. To be specific, let us consider (la); the case (lb) 
can be done by analogy. 

Let us quantize the directional volumes by means of Wi defined in (I2.2.8p . A 
standard procedure giveJl 

w := — - — ( O cosh (3 k A) + cosh (3 k A) O 

= -^(2cosh(&x)-^- + 6sinh(6a;)), (2.3.3) 

where a := ?>K^\fi and b := 3k, and where we have used the representation for the 
elementary observables defined by (12.3.11) . 

In what follows we solve the eigenvalue problem for the operator w and identify 
its domain of self-adjointness. 

Let us consider the invertible mapping L 2 {R,dx) 3 tp — > Uip —: f £ L 2 (I,dy) 
defined by 

U^(x):= ^ =:f(y), xeR, y G I := (0, tt/6). (2.3.4) 

sin 1 ' 1 (by) 



1 Subscripts of observables are dropped to simplify notation. 
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We have 



{ip\ip) = / ipip dx 

J — oo 

f ?(ln Itg 1 ^) |)^(ln |tg 1 ^(^)|) rf (ln |tg^(|) |)) 
f ^(l n |tgV &( ^)| ) ^(i ri | t gV &( ^ ) | ) _ ^ 



o 



2 1 & v 2 y sin(fa/) 

dy=(U^\U^}. (2.3.5) 



o 



sin 1 / 2 (6?/) sin 1 / 2 (fa/) 

Thus, the mapping ( 12. 3. 4ft is isometric and hence unitary. 

Now, let us see how the operator w transforms under the unitary map ( 12.3.4ft . 
The transformation consists of the change of an independent variable 

2 

x !->■ y := -arctan(e fex ), (2.3.6) 



which leads to 



la ( d \ d ab 

— — ( 2 cosh(fac)- — hfasinh(fac)) ^ —ia- — h i— cot(fa/), (2.3.7) 

Zi V (XJb / (-by ^ 

and re-scaling with respect to a dependent variable 

d . ctb . . . _i/2/7 \( . d ,ab , , \ , i/9 /T \ . d 
— ia- — h i — cotfcwj i — y sm ' (by) — ia- — h z — cot(by) sin '"(by) = —ia—. 
dy 2 V. */ 2 / dy 

(2.3.8) 

In the process of mapping 

w i-)> UwU' 1 = -ia— =: w, (2.3.9) 
ay 

we have used two identities: sin(fay) = 1/ cosh(fac) and sinh(fac) = — cot (fa/). 

Since to > (for O > 0), we assume that the spectrum of w consists of positive 
eigenvalues. To implement this assumption, we define w := and consider the 
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eigenvalue problem 

-a 2 ^-J v = v 2 f v , ye(0,n/b). (2.3.10) 

There are two independent solutions for each value of v 2 (where v e R), namely: 
sm(^y) and cos(^y). Removing this degeneracy leads to required positive eigenvalues 
of w. We achieve that in a standard way by requiring that the eigenvectors vanish 
at the boundaries, i.e, at y — and y = n/b. As the result we get the following 
spectrum 

f 1/ = Nsm(-y), v 2 = {nab) 2 , n = 0,l,2,... (2.3.11) 
d 

It should be noted that for n = 0, the eigenvector is a null state and thus the lowest 
eigenvalue is v 2 = (ab) 2 . Next, we define the Hilbert space to be the closure of the 
span of the eigenvectors ( 12.3. lip . The operator w 2 = —a 2 -^ is essentially self-adjoint 
on this span by the construction. Due to the spectral theorem [13] we may define an 
essentially self-adjoint operator w = \j—& 2 ^p. as follows 

wf u := vf u , v = ab, 2ab, 3ab, . . . (2.3.12) 



We have considered the case w > 0. The case w < does not require changing of 
the Hilbert space. The replacement w i-> — w leads to v \-> — v. 

Finally, we find that the inverse mapping from L 2 (I,dy) to L 2 (IR, dx) for the 
eigenvectors of w yields 



,v x , s ~ i , , , s sin ( 2 ^arctg(e 6a: )) 

sin -y = U(y) ^ fT 1 f u (y) := ^ (x = }) ■ 2.3.13 
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2.3.3 Kasner-like case 

In the case (2a), the conditions O1O2 + O1O3 + O2O3 > with d < 0, 2 > 0, 3 > 
are satisfied in the following domainsj for 

OiG(-di,0), 2 G(d 2 ,oo), 3 G(4,oo), (2.3.14) 

where d 2 > dx, and where ds = d\dij [d 2 — dx) so d 3 > d\. The full phase space sector 
of the Kasner-like evolution is defined as the union 

U {-d x , 0) x (d 2 , 00) x (d 3 , 00) (2.3.15) 

0<d 1 <d 2 

In the case of 2 and O3, the restrictions for domains (12.3. 14ft translate into the 
restrictions for the corresponding domains for the observables w 2 and u> 3 , due to 
(E2SD, and read 

w 2 G (£> 2 ,oo), W36 (^3,00), (2.3.16) 

where D 2 = fir/Ac^ and D3 = K/yAofa. Thus, quantization of the w 2 and W3 observables 
can be done by analogy to the Kasner-unlike case. The spectra of the operators w 2 and 
W3 are almost the same as the spectrum defined by ( 12.3. 121) with the only difference 
that now v > D 2 and v > D 3 , respectively^. 

The case of Wx requires special treatment. Let us redefine the elementary observ- 
ables corresponding to the 1-st direction as follows 

Q ^ : =-l T7TTT' CV=sinh(&A). (2.3.17) 

cosh (oAJ 

The transformation (12.3. 17ft is canonical, since {Qi,Q 2 } = 1, and invertible. The 
domains transform as follows 

O 1 e(-d 1 ,0), Ax eR — > fix G (0, =: (0,Di), fi 2 6l. (2.3.18) 



2 The case (2b) can be done by analogy. 

3 Spectra are insensitive to unitary transformations. 
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The observable V\ in terms of redefined observables reads 

vi = ^-Qi(l+0^), ^e(0,oo), (2.3.19) 

where ab/h = 127rG7A. To quantize observables fli and fl 2 we use the Schrodinger 
representation 

^2 -> n 2 /(x) := -ihd x f(x), -> Cltfix) := xf(x), f G L 2 (0, D^. (2.3.20) 

Let us find an explicit form for the operator tHOi + f^f^)) corresponding to 
(I2.3.19p . Since Qi > 0, the following classical equality holds 

Q^l = fij • 2 • fi^ fc - m • fi 2 • fi™, (2.3.21) 

where £ R. This may lead to many operator orderings at the quantum level. 
This issue is further discussed in the appendix. 

We propose the following mapping (we set h = 1) 

= -xd 2 xx -d x + rnkx' 1 , (2.3.22) 

which formally ensures the symmetricity of Q\Q?;. The second equality in ( 12. 3. 22ft 
may be verified via direct calculations. 

Now, we define the following unitary transformation W 

L 2 ([0,D 1 ],dx) 3 f(x) ^ Wf(x) := \[\f{^) G L 2 ([0, 2^}, dy). (2.3.23) 
One may verify that we have 

Wd x W* = -d v - \, Wd 2 xx W^ = \d 2 yy -\d v + \. (2.3.24) 

y y y y y 
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Thus, the operator W transforms (I2.3.22p into 



The eigenvalue problem for Oi + Vt\Vt\ reads 



- d w + — ( Amk - ~ J • (2.3.25) 



^ + ^( 4 ^-i) + xl $ = z/$ - (2 - 3 - 26) 



Now, we can see an advantage of the chosen ordering prescription (12. 3. 22ft . It enables 
finding a very simple form of the volume operator. Taking k = m = 1/4 turns (12.3.261) 
into 

- d lv + J " ") $ = °- ( 2 - 3 - 27 ) 

The problem is mathematically equivalent to the one dimensional harmonic oscillator 
in a "box" with an edge equal to 2\fD\. There are two independent solutions for a 
given v 

$u,x = N.e-y 2 " l F 1 ( - \v + ~ ~ ?Q , (2.3.28) 

= A^e"^ 4 i^i ( - ^ + y > y) , (2-3-29) 

where is a hypergeometric confluent function, $„ ! and $^2 are even and odd 
cylindrical functions, respectively. A standard condition for the symmetricity of the 
operator defining the eigenvalue problem (12.3.271) leads to the vanishing of the wave 



functions at the boundaries (as the box defines t 
The solution (after retrieving of h and ab) readc 



re entire size of the 1-st direction). 



$ = Nye~r h X F X ( - + -, -, 4 ] • (2-3.30) 
1 2ab 4'2'2fi' 



J We ignore the solution $^1 because it cannot vanish at y = 0. 
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The solution (I2.3.30P vanishes at y = as $ is an odd function. The requirement of 
vanishing at y = 1\[T)\ leads to the equation 

flu 3 3 2DA 



2abh 4' 2' h 

An explicit form of (I2.3.3ip reads 



n=0 V2/n \ / 

where (a)„ = a(a + 1) . . . (a + n — 1). It results from (12.3.321) that the eigenvalues 
must satisfy the condition: v > ab. 

2.3.4 Volume operator 

Classically we have 

V = \w 1 w 2 w 3 \ 1/3 . (2.3.33) 
One may verify that V/~ Poisson commute and dk commute, so we can take 

V 3 := V1V2V3 = \wiW2W 3 \. (2.3.34) 

The eigenfunctions of the operator wiw^Wz have the form 

^7 , Ai,a 2 ,a 3 ._ / 1 Al (x 1 )/ 2 A2 (^2)/3 A3 (^3), where fi Xi (xi) is an eigenvector of W{ with eigen- 
value Aj. The closure of the span of _F Al > A2 > A3 is a Hilbert space, in which V 3 is a 
self-adjoint operator (by construction). 

Due to the spectral theorem on self-adjoint operators [43J, we have 

V={V 3 ) l/3 — > yp^iM^s := □_pAi,A a ,A 3) (2.3.35) 

where 

□ : = lAxAaAal 1 / 3 . (2.3.36) 
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Kasner-unlike case 

In the Kasner-unlike case we use the formula (12.3. 12[) to get 

\3 = \riin 2 n 3 \ 1 ^ 3 ab, n 1; n 2 , n 3 G 1, 2, 3, . . . , (2.3.37) 

which shows that the spectrum of the volume operator does not have equally distant 
levels. The volume □ equal to zero is not in the spectrum. There exist a quantum of 
the volume which equals A := ah = l27rG'y\h, and which defines the lowest value in 
the spectrum. 

Kasner-like case 

The spectrum in this case reads 

□ := (J O dl42 , D dl42 := {X dl \ d2 X d3 \d 3 = d 1 d 2 /(d 2 - dt)}, (2.3.38) 

0<di<d 2 

where X dl is any value subject to the condition (I2.3.32[) . \ d , 2 > D 2 and A^ 3 > D 3 are 
given by ( 12.3. 12[) . The volume □ equal to zero is not in the spectrum. 

2.3.5 Evolution 

In this section we ignore the restrictions concerning the domains of Oi, O2 and O3, 
and we assume that the Hilbert space of the system is L 2 (R 3 , dxdydz). An inclusion 
of the restrictions would complicate the calculations without bringing any qualitative 
change into the picture of evolution. 

The generator of evolution determined in ( 12.2.131) may be formally quantized, due 
to ( 12.3.11) . as follows 

Hh>H = -h 2 (d y d z + d x d z + d x d y ). (2.3.39) 
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Since it is self-adjoint in L 2 (M. 3 ,dxdydz), a quantum evolution can be defined by an 
unitary operator 

U = e - ihT ( a y d z +dxd z +dxd y) ( r G R. (2.3.40) 
Let us study an evolution of the expectation value of the directional volume v\ 

(iplU^vxUlip) (2.3.41) 

Since V\ does not depend on y and z, we simplify our considerations by taking 

JJ X = e ~ ihT{dz+dy)dx . (2.3.42) 

If we are interested in the action of U\ on the functions f(x) G L 2 (R, dx), then 
the derivatives —i-^ and — i-^ occurring in U\ commute and, being self-adjoint, lead 
finally to real numbers. Let us call them k y and k z , respectively, and let us introduce 
the parameter k = k y + k z . Hence, U\ further simplifies and reads 

U x = e khTd *. (2.3.43) 

The action of U\ on f(x) reads 

UJ(x) = f(x + khr). (2.3.44) 

We recall that under the unitary mapping L 2 (R, dx) i-> L 2 (I, dy), defined by fl2.3.4p . 
the operator v\ becomes — ia-^ on L 2 (I, dy). Now, let us study an action of operator 
U\ on the functions <p(y) G L 2 (I,dy). Straightforward calculation leads to 

o, v , s v(l arctan(e 6a; )) „. . . 
L 2 (I, y) 3 V (y) ^ ^ U - JJ G L 2 (R, x), (2.3.45) 

cosh ' (ox) 

and we have 

arctan(e ba; )) u>(| arctan(e ba:+6fcriT )) 

^i— 1/2 = b 1/2 ~ (2.3.46) 

cosh ' (for) cosh ' (bx + bkhr) 
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The transformation U 1 gives 

^(§arctan(e 6a:+6fe?ir )) ^ ip{\ arctan(e fefc?lr tan(^))) _ ^ 

cosh 1 / 2 (6x + bkhr) R S i n (^)( tan (^)e^ + cot(^)e~ b ^) 



(2.3.47) 

where (p T=0 (y) = tp(y). Now, we observe that the symmetricity condition 

((p T (y)\vx<Pr(y)) = (vi<Pr(y)\<pM) (2.3.48) 

leads to 

Vt{\)Vt{\) - VMvM = 0. (2.3.49) 
We use the result (I2.3.47H to calculate the limits 

. . bkhr / \ i / \ bkhT , 7T > , , 

lira (p T (y)=e ? <p {0), km ip T {y) = e * <Po{-r), (2.3.50) 

y—tO y —+f b 

which turns (I2.3.49H into 



?o(£ )^{\)e- hkhT - v (0)Vo(0)e bkRr = 0. (2.3.51) 

It is clear that ( 12.3.51 j) can be satisfied Vr iff <^o(f ) = = <fio(0). States with such a 
property belong to the domain of w defined by (12.3.121) . 

In order to construct the "evolving states" that vanish at the boundaries, consider 
the basis vectors f n (y) = e l2bny . Then, f n (y) — f m (y) satisfy the condition ( I2.3.5ip . 
Making use of (I2.3.47H we get 



i - e bkhr tan(^)\ 2n / l + tan 2 (^ N 



/n ^' r) U + e 6M Uan(|)J \J + e bkhT tan 2 (^) ' ( 2 -3.52) 

where f n (y, T ) '■= fn,r(y)- Moreover we have 

- ia±f n (y,r) = -4(l + tan 2 (|)) /w (,,r) i + e2b J tan2( ^ x 



x ( " / , ^J,?^ + ^ne bhhT ) . (2.3.53) 



l + tan 2 (^ 
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Using the substitution x = tan(y) we get 



(/ml -ia-j-fn) 
dy 



x ' % — e bkhr x \ 2 ( n ~ m "> ( e -bkhr _ e bkhT\ x 







^ % I G bk ^ T X J (fy—bkhr _|_ ^bkhr j,2\2 (2.3.54) 



\i + e bkhT x) ( e -bkhr + e bkhr x 2\2 dx - 



Another substitution z = e T x leads to 



(fm\ ~ ia-^fn) 

dy 



-Me* -«•»)/ (^) JTTT^ dz ^ 

f°° ( { — z \ 2 ( n " m ) e bkHr _|_ e -bkhr z 2 

+4an / ( I ; — dz 

Jo V + z) (1 + ^ 2 ) 2 

Finally, we obtain 

d , I 77 — "\ 2 i (1 — 8n(n — m)) smh(bkhr), n^m 

(fm\-ia—fn) = { 4 (™- m ) 2 - 11 1 , \ ( 2 - 3 - 56 ) 

"2/ I zasinh(6/cftr) + 27macosh(6/i;frr), n = m. 

Now, let us introduce g n m(y, T ) '■= ^'^^^^^ so that ||y nm || = 1. One has 



(flWnl — ia-rQnm) = (w + m)abcosh(bkhT) = A coshibkhr). (2.3.57) 

dy 2 

The expectation value of the operator (T2.3.57]) . defining the volume operator, is 
similar to the classical form (I2.2.8P . The vectors g nm may be used in the construction 
of a basis of the space of states such that <^o(f ) — — ¥>o(0). 

2.3.6 Summary 

As in the FRW case, resolution of the singularity in the Bianchi I model is due to 
the loop modification of the Hamiltonian at the classical level. This modification is 
parameterized by a free continuous parameter A. 
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The spectrum of the volume operator, parameterized by A, is bounded from below 
and discrete. An evolution of the expectation value of the volume operator is similar to 
the classical case. We have presented the evolution of only a single directional volume 
operator. One may try to generalize this procedure to the total volume operator. In 
the case of the Kasner-like analyzes of dynamics are complicated. 

We introduced the so-called true Hamiltonian which proves an independence of 
the spectrum of the volume operator on the evolution. 



Conclusions 



Firstly, in the nonstandard LQC the results are obtained directly on the physical 
Hilbert space. Secondly, our nonstandard loop quantum cosmology, successfully ap- 
plied so far to the FRW and Bianchi I models, seems to be highly efficient. For 
example, only analytical calculations are needed to obtain the results. 

Turning the Big Bang into the Big Bounce in our method is due to the modification 
of the Hamiltonian at the classical level by making use of the loop geometry. The 
modification is parameterized by a continuous parameter A, which value is not known. 

In both considered models the spectrum of the volume operators, parameterized 
by A, are bounded from below and discrete. An evolution of the expectation values 
of the volume operators are similar to the classical cases. 

In the Bianchi I model, the phase space of the system is divided into the two 
distinct regions: the Kasner-like and the Kasner-unlike. Domains, spectra and eigen- 
vectors of self-adjoint directional volumes, and total volume operators were identified 
in the Kasner-unlike case. The peculiarity of the Kasner-like case was identified due to 
complicated boundary of the phase space region. We propose to overcome this prob- 
lem by dividing this region further into smaller regions, but with simpler boundaries. 
Given a small subregion for the Kasner-like case, we propose a canonical redefinition 
of phase space coordinates in such a way, that we can arrive at relatively simple form 
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of volume operator and at the same time can simply encode the boundary of the re- 
gion into the Schrodinger representation. Then, from a number of different operator 
orderings the simplest one was chosen. Domain, spectrum and eigenvectors of the 
volume operator was founded. The spectrum is given in an implicit form in terms of 
special functions. 

Discreteness of space at the quantum level may lead to a foamy structure of 
spacetime at the semi-classical level. The discreteness is also specific to the FRW 
case [38J. The difference is that in the Bianchi I case the variety of possible quanta 
of a volume is much richer. On the other hand, the Bianchi type cosmology seems 
to be more realistic than the FRW case, near the cosmological singularity. Thus, an 
expected foamy structure of space may better fit cosmological data. Various forms of 
discreteness of space may underly many approaches in fundamental physics. So its 
examination may be valuable. 

As we know A is a free parameter. Without specific choice of A, the Big Bounce 
may occur at any low or high density. The former case (big A) contradicts the 
data of observational cosmology (there was no Big Bounce in the near past) and 
leads to weakly controlled modification of the expression for the curvature F^ b , i.e. 
gravitational part of the Hamiltonian (see Appendix A). On the other hand the latter 
case (small A) gives much better approximation for the classical Hamiltonian (see 
Appendix A), but may easily lead to densities much higher than the Planck scale 
density, where the classical formalism is believed to be inadequate. Finding specific 
value of the parameter A, i.e. the energy scale specific to the Big Bounce is an open 
problem. It may happen, that the value of the parameter A cannot be determined, 
for some reason, theoretically. The story may turn out to be similar to the case 
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of the short-range repulsive part of the potential of the nucleon-nucleon interaction 
introduced to explain the scattering data [54] and the nuclear matter saturation of 
energy [55] . In such a case A will become a phenomenological variable parameterizing 
our ignorance of microscopic properties of the Universe. Fortunately, there is a rapidly 
growing number of data coming from observational cosmology that may be useful in 
this context. The cosmic projects for the detection of gamma ray bursts may reveal 
that the velocity of cosmic photons depend on their wave lengths, which may be 
ascribed to the foamy nature of spacetime [SB1 E3, EE]- Such dependence is weak, 
but may sum up to give a measurable effect in the case of photons travelling over 
cosmological distances across the Universe [2H]- Presently, available data suggest 
that such dispersion effects do not occur up to the energy scale 5 x 10 17 GeV [60] 
so such effects may be present, but at higher energies. Another way to determine 
the phenomenological value of the parameter A is the detection of the primordial 
gravitational waves created at the Big Bounce [HI (62J [631 El]. 

In our method an evolution parameter (p does not belong to the physical phase 
space, contrary to the standard LQC Thus, it stays classical during the quantization 
process as well. At quantum level of the Bianchi I model, for the first time in our 
method, the so-called true Hamiltonian was introduced. It generates a flow in the 
family of volume quantities, enumerated by an evolution parameter. Having the true 
Hamiltonian, we could introduce an unitary operator with the evolution parameter 
t el 

It is clear that the next step, in the road to understand an early Universe, is the 
nonstandard LQC quantization of the Bianchi II cosmological model. 



Appendix A 
Holonomy corrections 



The curvature of SU{2) connection F k b = d a A b — dbA^ + e^A^Al, entering the expres- 
sion (11.2.2P for the gravitational part of the Hamiltonian, can be expressed in terms 
of holonomies. Using the mean- value and Stokes' theorems we have 

4- / T k F* d dx c Adx d n-^(vexp(<f r k A k c dx c )-l), (A.0.1) 

S ab J<J S ab V J da J 

where da is the boundary of a small surface a with center at x, and where s a ab := 
J dx a A dx b . The expression for F k b is exact but in the limit when we shrink the area 
enclosed by the loop da to zero. If we choose da in the form of the square with 
sides length A, the expression for a small value of A = /io has the form [65] 

F "'<"»> = a™ { - 2 Tr (fi^) r " °< °< + < A -°- 2 > 

and we have 

F* = lim F* (/Iq). (A.0.3) 



In the standard LQC the 0(X A ) holonomy corrections are ignored (see, e.g. 
). It was found in [521 EE] that including higher order corrections leads to new 
curvature singularities different from the initial singularity and increases an ambiguity 

60 



61 



problem of loop cosmology. However, the holonomy corrections do not change the 
result that the Big Bounce is a consequence of the loopy nature of geometry [67] . 

Taking only the first term of (lA.0.2j) leads to the simplest modification of gravity, 
but may be insufficient for the description of the inflationary phase. The choice of /xq 
based on the expectation that the Big Bounce should occur at the Planck scale [22] has 
little justification [SB]. The significance of Planck's scale for quantum gravity seems to 
be rather a belief than proved result (see, e.g. [68J). Heuristic reasoning playing game 
at the same time with Heisenberg's uncertainty principle, Schwarzschild's radius and 
process of measurement cannot replace a proof (see, e.g. |69|). 



Appendix B 

Justification for using observables 
to classify phase space 



B.l Nonregularized case 

The nonregularized case means that we do not modify general relativity by loop 
geometry (loops of finite length). For this case one haql 

0, = —Vipi, (B.l.l) 
^7 



where 



It is clear that 



which leads to 



Vi := (a^L^ 2 , ft := °f- (B.1.2) 



Oi = —aidjakLiLjLk— =: — Hi, (B.l. 3) 



H k = ^O k (B.1.4) 

where H k '■— a k ja k is a directional Hubble parameter in the /c-th direction. In what 
follows we use H k to describe contraction or expansion of spacetime. It results from 



1 obtained from (|2.2.1j) in the limit A -> 
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(1B.1.4j) that classical dynamics of our system in the nonregularized case can be de- 



scribed in terms of the observables O k . 

B.2 Regularized case 

In this case we modify gravity by loop geometry, due to (12.2. 6p . we have 



Vi = K-f\\Oi\ cosh (3K((Oj + O k )<p- A*)) (B.2.1) 

so we get 

— = 3k (Oj + O k ) tanh (3k((Oj- + O fc ) <p - A;)) (B.2.2) 
where (p = ~ ^. On the other hand we have 

V -i = h H] + H k ). (B.2.3) 

Comparing (1B.2.2j) and (IB.2.30 we obtain 

Hj + H k = ^p (Oj + O fc ) tanh (3/t((O i + O fe ) ^ - A,)). (B.2.4) 

Denoting the moment of occurring a bounce by <£>f , defined by 

A i = {O j + O k )i P f, (B.2.5) 

we obtain the formula 

Hj + H k = £ (Oj + O k ) tanh (3k (O, + O fc ) (<p - <pf)). (B.2.6) 

For the purpose of interpretation it is convenient to rewrite (IB. 2. 61) in different 
form. Since cosh (x) 2 — sinh(x) 2 = 1, we have 
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In the case of expanding universe, ip > iff , we have the expression 



tanh(3K(O i + O fe ) {</>-$)) = Jl 



It results from (IB.2.ip and flB.2.5|> that 



cosh 2 (3k {Oj + O k ) (if - iff)) ' 



(B.2.8) 



n^XOi cosh (3k (Oj + Ofc) (if — iff)) ■ 



Since vf = fi/yAO, (directional volume at the i-th bounce), we have 



(B.2.9) 



cosh (3k (Oj + O fc ) (if - iff)) 

which means that 

tanh (3k (Oj + O fc ) (if - iff)) = 
Inserting (|B.2.11[) into flB.2.6[) gives finally 



Vi 



(B.2.10) 



Vi 



(B.2.11) 



(B.2.12) 



It results from flB.2.12D that for very large Vi, comparing to vf, we can apply the 
approximation 



Hj + Hu 



K 

V 



(Oj + A 



(B.2.13) 



which finally leads to 

H k = ^O k . (B.2.14) 

Comparing (IB.1.4[) and (IB. 2. 141) one can see that also for the regularized case the 
structure of the classical phase space may be described in terms of observables, 
but only when we consider the Universe far away from the Big Bounce (vi 3> vf). 
However, for our analyzes that limitation does not matter. 



Appendix C 

Non-uniqness in quantization of 
the Bianchi I 

It is clear that quantum cosmology calculations are plagued by quantization ambigu- 
ities. For example, there exists a huge freedom in ordering of elementary operators 
defining compound observables, which may lead to different quantum operators. Clas- 
sical commutativity of variables does not extend to corresponding quantum operators. 
Other ambiguities are discussed below. Such ambiguities can be largely reduced when 
some quantum data from cosmological observations become available. Confrontation 
of theoretical predictions against these data would enable finding realistic quantum 
cosmology models. 

C.l Unitarily non-equivalent volume operators 

In both Kasner-like and Kasner-unlike cases, we have reduced the Hilbert space by 
removing the double degeneracy of eigenvalues for the volume operators (see the dis- 
cussion after equations (12.3. 10p and (12.3.27]) ) . We have used the "natural" condition 
that the wave function should vanish at the boundaries of an interval. However, there 
are also other mathematically well-defined choices for the boundary conditions. We 
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will demonstrate this non-uniqueness for the Kasner-unlike case. Similar reasoning 
applies to another case. 

Let us begin with the equation (12.3. 10p 

-a 2 ^/ = ^ 2 /, ve{0,*/b), (C.l.l) 

which has the solution 

/„ = m sin(-y) + N 2 cos(-y), N 1} N 2 e C, (C.1.2) 
a a 

for each value of v G M + ( z/ i-> — 1/ does not produce any new space of solutions). 



Our task is the determination of self-adjointness of w := y — fl2 J^ and removing the 
double degeneracy of eigenvalues. The symmetricity condition reads 



/» Tt/b n / u p 

/ //" = //' " f'f + / /"/• 

7/ o o Ji 

We can set: 

• /(0) = /(tt/6) = ^ /„ = sin(^), v = ab, 2ab, 3ab, . . . 

• f(0) = /'(tt/6) = f v = cos(^), v = 0, 06, 2a6, 3a6, . . . 

• /(0) = /'(tt/6) = =}► /„ = sin(^), i> = iafe, fafc, f ab, . . . 

where ab = 127ihG^f\. All these choices are non-equivalent, since they lead to different 
spectra. 

C.2 Standard quantization 

Let us change the coordinates of the Kasner-like sector phase space Q 2 ), defined 
by (12.3. 17p . into a new canonical pair as follows 

X := and P := ^2^2^, (C.2.1) 



ir/b 



(C.1.3) 
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where 

(X, P) G (0, yfidk/b) x R, {X, P} = 1. (C.2.2) 
In the new variables the volume (12.3. 19j) reads 

' Vl = -P 2 + ix 2 . (C.2.3) 



4ttG , 7A 2 2 

Thus, in these variables the volume has a form of the Hamiltonian of the harmonic 
oscillator in a "box" (0, yj2d\/b). 

In the Schrodinger representation, i.e. X := x and P := —ihd x , a standard 
quantization yields 

1 v = -%d\ x + \x\ (C.2.4) 



47rG 7 A 2 3535 2 

which corresponds to the "nonstandard" quantization f|2.3.26[) with the parameters 
m = k = 1/4 and y = \/2x (with h— 1). 

Thus, we can see that the prescription defined by (I2.3.2ip and (I2.3.22j) includes 
not only a standard prescription, but many others. As an illustration only one, 
corresponding to the well known harmonic oscillator, has been completed. 
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